Charge transport in graphene with resonant scatterers 



M. Titov,i'2 p. M. Ostrovsky,3>4 i. y. Gornyi,3>5.2 Schuessler,^ and A. D. Mirlin^' '^^ 2 

' School of Engineering & Physical Sciences, Heriot-Watt University, Edinburgh EHI4 4^S, UK 
^ DFG Center for Functional Nano structures, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany 
Institut fiir Nanotechnologie, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany 
L. D. Landau Institute for Theoretical Physics RAS, 119334 Moscow, Russia 
' A. F. loffe Physico- Technical Institute, 194021 St. Petersburg, Russia. 
^ Institut fiir Theorie der kondensierten Materie, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany 
^ Petersburg Nuclear Physics Institute, 188300 St. Petersburg, Russia. 

The full counting statistics for the charge transport through an undoped graphene sheet in the 
presence of strong potential impurities is studied. Treating the scattering off the impurity in the s- 
wave approximation, we calculate the impurity correction to the cumulant generating function. This 
correction is universal provided the impurity strength is tuned to a resonant value. In particular, 
the conductance of the sample acquires a correction of 16e^/(7r^/i) per resonant impurity. 

PACS numbers: 73.63.-b, 73.22.-f 



Since the discovery of graphene [H its transport prop- 
erties have become a subject of intense studies The 
most remarkable effects arise when the chemical poten- 
tial is tuned into a close vicinity of the Dirac point. In 
particular, a short and wide sample of clean graphene 
exhibits a pseudo-diffusive charge transport Q , with the 
counting statistics equivalent to that of a diffusive wire 
[^,^] . This equivalence has been confirmed in recent mea- 
surements of conductance and noise in ballistic graphene 
flakes In contrast to conventional metals, ballis- 

tic graphene near the Dirac point conducts better when 
potential impurities are added (9l-[lTl| . Quantum interfer- 
ence in disordered graphene is also highly peculiar due to 
Dirac nature of carriers. In particular, in the absence of 
intervalley scattering, the minimal conductivity /h 
is "protected" from quantum localization 

Strong impurities creating resonances near the Dirac 
point ("midgap states") are one of the most plausible 
mechanisms limiting the electron mobility and can be 
used exploited for functionalization of graphene. As was 
shown in Refs. [l3|,|lj], such scatterers provide the con- 



centration dependence of the conductivity, cr cx nln n, 
which is consistent with most of experimental observa- 
tions. Possible realizations of such scatterers are va- 
cancies, adsorbed atoms, molecules, or impurity clusters 
[isl [l6|. In particular, important example is hydrogen 
atoms that can be controllably added to the graphene 
sample [l7'|. Resonant scattering may also be intention- 
ally induced by metallic islands deposited on graphene 
surface [3] . In this case the strength of impurity can be 
controlled by a local external gate. 

In this paper we study the effect of strong impurities 
on electron transport in graphene. We consider the bal- 
listic transport regime in which the sample size is smaller 
than the electron mean free path (low impurity concen- 
tration) . This model was experimcntallyimplcmcntcd in 
the samples on the Si02 substrate 0,0, [3]. The ballistic 
transport is particularly relevant for suspended samples 




FIG. 1: (Color online) Ballistic graphene setup with various 
strong scatterers. Vacancies as well as atomic or molecular 
impurities can create midgap states 15]. Metallic islands sup- 
port quasibound states that can be tuned to the resonance. 



where much higher mobilities have been achieved 2^, 21 1 . 

We evaluate analytically the effect of strong impuri- 
ties on the full counting statistics (FCS) in a rectangular 
sample of length L and width W (Fig. [Ij. In particular, 
for the short-and-wide setup, ly L, we find the uni- 
versal corrections to the conductance G and to the shot 
noise S (divided by 2eV, where V is the bias voltage). 
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per resonant potential scatterer of the round shape. 

We consider an effective model described by the Dirac 
Hamiltonian, H — —ihv(T'V + U{r), where cr = {ax,o-y) is 
the vector of Pauli matrices, v is the velocity, and U{r) 
is the impurity potential. Below, we set hv — 1. Two 
metallic leads at cc < and x > L are defined by adding 
a large chemical potential to the graphene Hamiltonian. 
Inside the sample, i.e. for < x < L, the chemical poten- 
tial is set to zero (Dirac point). The potential U{r) rep- 
resents a collection of randomly distributed sharp scat- 
terers of arbitrary strength. 

We will study FCS of ballistic electron transport 



2 



through the sample using two complementary approaches 
— scattering theory and Green function formalism p^ . 
The first approach works for any aspect ratio of the sam- 
ple while impurities are treated as point-like. The second 
method is particularly suitable ior W ^ L and allows for 
arbitrary impurity profile. 

We begin with the scattering approach for a single im- 
purity We discretize y coordinate, y = Wv/N, where 
I' = 1, 2, . . . , is an integer, N = W/ na^ 1, and a <^ L 
plays the role of the impurity size. The impurity poten- 
tial in this model is defined as 



U(r) = a d{x — xq) 6^, 



(2) 



where Xq and yo = Wvq/N specify the impurity posi- 
tion. Generically, a is a 4 x 4 matrix in sublattice and 
valley spaces. Below we concentrate on the case of scalar 
potential with a being a number. 

In order to study transport properties, we perform a 
standard unitary rotation of the Hamiltonian, H i— 
C'^HC with C — {(Jx+<^z)/V^- After this transformation 
the upper (lower) element of the spinor wave function 
represents right- (left-) propagating mode in the leads. 
We then perform the discrete Fourier transform with re- 
spect to V and arrive at the Dirac equation 



d^{x) 
dx 



a^q - ia^U{x) *(x). 



where g is a diagonal matrix with entries qn = 
2nn/W being the discrete transverse momenta and n = 
—N/2,..., N/2 (for definiteness, we assume periodic 
boundary conditions). The impurity potential is repre- 
sented by the operator U{x) — aS{x ~ xo)\^){^\ that 
projects onto the state with the wave function <&(<7n) = 
^-i/2gjij„yo ^ Using this separable form of the impurity 
potential U, we can explicitly calculate the transfer ma- 
trix T that relates the wave-function amplitudes in the 
opposite leads, ^'(L) = T5'(0). The result is given by 

g(T^(j(L-a;o)gjcr^ 

Inverting the element Tn (where indices refer to the a 
space), we obtain an exact expression for the transmis- 
sion amplitude from the nth channel in the left lead to 
the mth channel in the right lead (see Appendix E| . 



27„,„[z-f iiVcot(a/2)]-i 



(3) 



coshg„L cosh(7„L coshg,„L 

with 7„,„ = g<q'^-<i^)yo cosh[(7„xo - qm{L - xq)] and 

z = ^ cosh[(j„(i - 2xo)]/ cosh((7„L). (4) 

n 

The first term in Eq. ([3]) describes the clean system 
{a — 0). The second term represents the effect of the 
impurity which is particularly important in the vicinity 




W/L 

FIG. 2: (Color online) Contribution to the conductance G 
(solid) and to the shot noise S (dashed) from a resonant im- 
purity placed at xo = L/2, as a function of the aspect ratio 
W/L. Dotted lines show asymptotic values from Eq. ([T]). In- 
set: conductance and noise of a clean sample, cf. Ref. [5|. 



of resonant values a„ = 7r(2n-f 1). The quantity z defined 
by Eq. Q is nW times the local density of states of a 
clean system at the position of impurity. 

The statistics of the charge transport is described by 
the cumulant generating function T{x) as c„ = 
limp^_j,o where ci —G/g^) is the dimensionless 

conductance, ci = S/2eVgQ is the dimensionless noise, 
etc. For our purposes it is more convenient to use the 
variable defined via x = 2 In | cos((/i/2)|, which yields 



F{(t)) ^ Indet [I - W sin2((/)/2)] . 



(5) 



Using Eq. ([3]) , we calculate conductance and noise of the 
sample with a resonant impurity placed in the center, see 
Fig-Hi The general expression for the generating function 
applicable for any strength and position of the impurity 
is given in Appendix VK\ 

In the limit W ^ L the summations over transversal 
modes can be replaced by momentum integrals. For the 
clean system this yields J-q = — M^0^/47rL. The impurity 
correction to J- takes the form 



(5 J" = 2 In 



I - 



1 



(4L)2 V7r2 sm^{TTXo/L) 



(6) 



where £ — 27ra tan(Q;/2) (as we show below, £ is the scat- 
tering length at low energies). In particular, the correc- 
tion to the conductance Go — goW/irL reads 



SG = {8go/TT^) [sin-2(7rxo/i) + (il f/l^ 



(7) 



This completes our analysis of the discretized model © . 

Let us now turn to the microscopic description of 
the impurity within the continuous Dirac model, assum- 
ing that the impurity potential is rotational symmetric, 
U{r) = M(|r — TqI), where rg is the impurity position. To 
calculate FCS we use the matrix Green function approach 
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[6|, |22|. The Green function in the retarded-advanced 
(RA) space satisfies the equation 



L) ii{x)-H-iQ 



G(r,r') = 5{v-y') 



where C = sin (0/2) is the counting field and /i is the 
chemical potential which is zero inside the sample and 
infinite in the leads. An explicit solution to the above 
equation is given in Appendix [B] 

The generating function can be expressed through G 
as = TrlnG"^, where Tr is the full operator trace. 

Expanding T{4>) in U{t), we obtain the impurity correc- 
tion to the generating function as a series 



OO -. 



(8) 



where Go is the Green function of the clean system. Eval- 
uation of Eq. ([HI essentially simplifies if the impurity size 
a is small compared to L. In this limit the arguments of 
the Green function are close to the position of the impu- 
rity, ro . Therefore it is useful to split Gq into the singular 



and regular parts, Gq 



Greg, where 



g(r,r') = -(*/2^V.(r-r')/|r-rr 

is the zero-energy Green function of the free Dirac 
fermion. The regular part of the Green function, Greg, 
can be taken with equal arguments r = r' = Vq. Ex- 
plicit expression for Greg(r) is given in Appendix |B] We 
characterize the impurity by its T-matrix 

OO „ 

T^Y. ^'^1 ■ • ■ c?'rp[/(ri)5(ri, ra) . . . U{vp). 

p=i-' 

Rearranging singular and regular parts in Eq. ([8]) yields 

(5J- = Trln[l-TG,eg(ro)]. 

Note that this result holds for any T matrix, includ- 
ing possible valley mixing. In the case of rotational in- 
variant impurity, the T matrix reduces to the scattering 
length ^ in the s-channel. Since TrGiog(ro) = and 
G?eg(ro) = [0V^'-sin-2(7ra;o/L)]/(4i)2 (Appendix |B| , 
the impurity correction to the generating function repro- 
duces the result ([6]). This establishes a relation between 
the phenomenological parameter a in Eq. ^ and the 
actual profile of the impurity potential. 

It is convenient to derive the scattering length t from 
the solution of the Dirac equation in an infinite system 
[2^ 1 at a low energy k. Specifically, i is related to the 
scattering cross-section A as A = M'^ /2. As an example, 
let us calculate I for an impurity with the rectangular 
profile u{r) = uo9{a — r). Matching the solutions of the 
free Dirac equation at r < a and r > a, we obtain 



£ = 2Tra JiiuQa)/ Jo{uQa), 



(9) 




FIG. 3: (Color online) Correction to the conductance from 
a circular impurity at xq = L/2 as a function of uoa. Solid 
line: analytic result, Eqs. (O, (O; dotted line: numerical 
simulations [2^ with the parameters a/L = 0.2, W/L = 6. 

where J„ stands for the Bessel function. The scatter- 
ing length diverges when the impurity potential develops 
a quasibound state (one component of the spinor wave 
function is localized) at zero energy. This implies reso- 
nant scattering, which is a generic feature of any strong 
impurity independent of its profile. In the quasiclassical 
limit, uqu ^ 1, one finds £ — 27ratan(uoa — 7r/4), which 
corresponds to a = 2uoa — 7r/2 in Eq. 

In Fig. Owe compare the results (O, © to the conduc- 
tance calculated numerically in Ref. [25] for the case of 
the disk impurity. In addition to the s-wave resonances, 
that are perfectly described by our theory, there exist 
higher resonances that can be included into our consid- 
eration in a similar manner, see Appendix [Cl The posi- 
tion of the resonances are given by the zeroes of Jm(uoa), 
where the corresponding index m = 0, 1, 2, ... is specified 
in Fig. m The resonance widths scale as (a/L)^'"+^. 

Now we turn to the case of disordered graphcne with 
a small concentration of impurities randomly distributed 
over the sample. The average generating function can 
be found to the linear order in nimp (see Fig. S^) by the 
integration of Eq. ^ over the sample area. 



T = — 



mn.p In 1 



1 



16L2 



where Ni^p — rii^pWL is the total number of impurities. 
Resonant impurities {£ = 00) contribute on average the 
universal correction ([T} to the conductance and noise, 
yielding the conductivity a = GL/W and Fano factor 



1 



12 



3 TT 

Away from a resonance, a ^ £ L, we find 



2 ^imp-^ 



■nh 



27r 



F 



(10) 



To establish the limits of validity of these ballistic re- 
sults, we estimate the contribution of the order nf^p in 
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comes into play. The symmetry-breaking length Lc is 



given by 



FIG. 4: (Color online) Conductance vs. L for resonant im- 
purities. Dotted lines: (a) one- and (b) two-impurity contri- 
butions. Inset: corresponding typical diagrams. Dasiied line 
illustrates schematically the crossover to diffusive regime. 

the virial expansion. The diagrams involving two reso- 
nant impurities, Fig. |4}d, yield the free-energy correction 

5(2)j-^ lTrln[l-G„i(ri)Go(ri,r2)G-i(r2)Go(r2,ri)], 

which is a contribution of atypical pairs of scatterers sep- 
arated by a distance < L. The corresponding corrections 
to the conductance and Fano factor are (Appendix |D]) 



J^^V = -1.032— n?L^ 
h ^ 



^ 0.8623 nf,„pL^ 



Two closely located strong impurities detune each other 
from the resonance, hence the negative correction to the 
conductance. For off-resonance impurities, £ L, the 
multiple scattering correction to conductance can be de- 
rived from the first-order result Eq. ([TUl) supplied by 
the logarithmic renormalization of the effective impurity 



strength ao 



4e2 



J^/2tt according to Refs. 10, 13], 



1 + 



2[7r - m^pi^ log(L/£)] 



When Sa ^ /h (which corresponds to mean free path 
^ L), the system enters the diffusive regime. Fig. 01 This 
happens at n\-^p min{L^, P\ log(L/£)|} ~ 1. Remarkably, 
for resonant impurities this implies just one impurity per 
square Lx L. The scaling of conductivity in the diffusive 
regime is determined by the symmetries of impurities in 
sublattice and valley space [12|. In particular, for vacan- 
cies the quantum interference is suppressed in view of 
chiral symmetry [l^l . This should yield the universal (as 
long as iVimp is small compared to the number of carbon 
atoms) minimal conductivity '-^ e^//i. 

The situation with strong potential scatterers (e.g., ad- 
sorbants inducing midgap states) is similar for interme- 
diate L where the chiral symmetry is approximately pre- 
served, Fig.m At larger scales, L > L^, violation of this 
symmetry becomes important, and quantum interference 



-1/2 
''imp 



minjnimp^^, l/n-impO^}: the chiral 



symmetry is broken due to detuning from resonance as 
well as due to higher (non-resonant) scattering channels. 
For impurities that do not mix the valleys (i.e. smooth on 
the scale of the lattice constant) this implies antilocaliza- 
tion (symplectic symmetry class) leading to logarithmic 
increase of conductivity with L. The valley mixing re- 
stores conventional (orthogonal) 2D localization. 

In conclusion, we have computed FCS for charge trans- 
port through an undoped graphene sheet in the presence 
of strong potential impurities. The impurity correction 
to FCS, Eq. (I6|), is determined by the position of impurity 
and the low-energy scattering length. The latter diverges 
when the impurity potential has a quasibound state at 
zero energy. At such resonant conditions the impurity 
correction becomes universal, Eq. ([1}. Our results are 
fully supported by numerical simulations. Fig. [31 
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Appendix A: Transfer matrix approach 

The full counting statistics of the charge transport is encoded in the cumulant generating function given by the 
formula 

= Indet [1 - W sin2(0/2)] , (Al) 

where t is an x iV matrix of transmission amplitudes, which are defined with respect to the scattering states in the 
leads. The matrix t is obtained by inverting the block 7ii of the transfer matrix T. 

For the case of a single scalar impurity, the 2N x transfer matrix is given by the product of three matrix 
exponents 

7" = gO";cg(i-a:o)giaCT^ | *) {$| ^cr^ gajo (A2) 



where the coordinate xq specifies the position of the impurity in x direction. The expression (jA2|) relays upon the 
sharp boundary conditions at the metal-graphene interfaces: a; = and x = L. For periodic boundary conditions in 
y, the momentum quantization gives = 7r(2n — N — 1)/W, where n = 1,2, ... ,A'^ (for simplicity, the number of 
channels, is taken to be an odd number), and the matrix q is defined as q = diag(gi, ^2, 13, ■ ■ ■ , Qn)- It is already 
assumed in the construction of Eq. (|A2[) that the impurity is point-like in x direction. Furthermore, if the impurity is 
point-like in y direction the vector |<i>) has the elements ^{qn) — A^~i/^e**"^° , where yo specifies the impurity position 
in y direction. The impurity strength in this phenomenological model is characterized by the dimensionless parameter 
a. 

In the limit of point-like impurity one can use the projective property of the operator |<i>)(tf>| in order to write 

e"""' = 1 - (1 - e"^"') |$)($|. (A3) 

Then the expression (IA2[) for the transfer matrix can be recast in the following form: 

where we have introduced the abbreviation D = 1 — exp(zao'2). It is also convenient to take advantage of 2 x A 
matrix 

^ _ 1 /coshgix coshg2a; coshqaa; ... coshg7va;\ ^^y^ (A5) 
^ ^/N \s'mhqix sinhg2a; sinhgsa; ... smhqNxJ ' 

which helps rewriting the matrix of transmission amplitudes, t ~ (7ii)^^, as 

* " coshqL-Wl_^^DW,, ~ coshqL + cosh<7i^^--«i5-i _ w,,^^Wl~^^^"" coshqL ^^^^ 

Thus the inversion of A^ x A^ block of the transfer matrix is reduced to the inversion of 2 x 2 matrix in the last 
expression. Indeed, the straightforward calculation yields 



cosh qL^ ^-^0 ' 2 N^V' cosh<z„L ' 



Substituting this result into Eq. (IA6[) we obtain 



11+1 2 2A'cr;, 

W; KW^ K = = 

coshgL coshgL "coshgL' 1 + (z/Af)^^ - 2/?-! z - i A cot(a/2) ' 



This expression is equivalent to Eq. ([3]). Using Eq. (jASp one can write the matrix product iH as 

= —T-^ + W^^W) -^-rr, (A9) 
coshgL V / coshgL 

where is a rectangular matrix of the dimension 4 x A" and 4x4 matrix S is responsible for the impurity correction 
to the transmission probabilities. 



W=\^^^ 1 I , S = I ° cosh^ qL ^-^^ I . (AlO) 

-A' 



Wl- 



' cosh qL / 
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Using the decomposition (jA9p we can essentially simplify the calculation of the determinant in Eq. (IA1[) for the full 
counting statistics. First, we factorize the determinant into the product of two determinants: one for the full counting 
statistics of a clean system and the other one for the impurity correction, 



det 1 



cosh^ qh ) 



det 1 



sin2(0/2) 



cosh qL — sin (0/2) 



(All) 



The first determinant of the diagonal matrix is evaluated straightforwardly. In the second one we perform a cyclic 
permutation of matrices and reduce it from N x N down to 4 x 4. This matrix is diagonal in a space and hence 
the determinant further reduces to the product of 2 x 2 determinants. After some algebraic manipulations we find 
out that these two latter determinants are identical. This allows us to represent the generating function as a sum 
T = To + 5T , where describes the clean system and 6T provides the impurity contribution to the full counting 
statistics. 



1 - 



sin2(0/2) 



cosh qnL 



(5 J" = 2 In 



1 



P(0) - Rl sin^ (j) 



cot2(a/2) 



Here we have introduced the following notations: 



Rx — ^ 



cosh(2g„a;) 



2 ^ cosh\qnL) - sm\(f>/2) 



PW = Rl-,, + 2Rl-,„R,o cos( 



(A12) 



(A13) 



The generating function T = To+SJ-' given by Eq. (|A12p is a general solution of the single- impurity problem applicable 
for any strength and position of the impurity and arbitrary aspect ratio of the sample. 

In the limit W ^ L the summations over transversal modes can be replaced by momentum integrals, so that 



W sm{(j)x/L) 
Rx — ■ 



2L sm{Trx/L)sm(j)' 
As a result, one obtains Eq. ([B]). 



Rf) — 



W(t> 



W 



2'kL sin ( 



2L sin(7ra;o/i) ' 



(A14) 



Appendix B: Matrix Green function 



The full counting statistics of the electron transport is conveniently expressed in terms of the matrix Green function 
[2^ in the external counting field C = sin(0/2). For the clean graphene sample this Green function satisfies the 
following equation in the retarded-advanced space: 



li{x) — crp + iO ~(JxC,5{x) 
^<^xC5{x — L) /i(x) — crp — iQ 



Go(r,r') = <5(r-r'), 



^l{x) 



0, < X < L, 

+00, X < or a; > L. 



(Bl) 



Since the operator in the left-hand side of the above equation commutes with the y component of the momentum, 
we will first calculate the Green function in the mixed coordinate-momentum representation, Gjy{x,x'). Inside the 
sample this function satisfies 



d 



Gp{x, x') — S{x — x'). 
We will look for a general solution of this equation in the form 



Gp{x,x') = e'^^f(=^-i/2)Me'^^''(^'~^/2)^ 



M = 



Af<, X < x' , 
My, X > x'. 



(B2) 



(B3) 



The chemical potential profile together with the infinitesimal terms ±iO in Eq. (|B1|) defines the boundary conditions 
for the Green function. The counting field C can also be incorporated into the boundary conditions. In terms of 
we thus obtain 



1-1 



g-^3pL/2^,^^ = 0, 



1 -10 

-iC -< 1 1 



g<T,pL/2^,^ = 0. 



(B4) 
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Delta function in the right-hand side of Eq. (jB2[) yields a jump of the Green function at x = x' which provides the 
relation 



The matrices Af^, and hence the Green function, are completely determined by Eqs. (jB4| IB5I 

/ coshpL ,r.-p^ 



(B5) 



M< = ^ 

> 2{cosh^ pL-C) 



^2 sinh 2pL 



\ 



cosh pL 

K 



■ coshpL 



2 sinh 2pL 



2 ^ sinh 2pL 



± 



(B6) 



cosh pi y 



Fourier transform in p yields the Green function in the full coordinate representation. To facilitate further calcula- 
tions, we decompose this Green function into the following product of matrices: 



Goix,x';y) = -V{x)A 



1 1 cosh — smh — 
2L 2L 



I 



sinh 



<\>y 



A = 



2L 

(Jz 
1 



-i cosh 



1 



1 



sin ^ {x + x' + iy) sin -^{x — x' + iy) 
1 1 



2L> RA Vsin^(.T-x'-iy) sin ^ (a; + x' - i y) / ^ 
V{x) = 



kV-\x% (B7) 



<b(L — a;) 6{L — a;) ' 
sm — cos ■ 



I cos ■ 



2L 

(f)X 

2L 



I sm ■ 



2L 

(j)X 

2L 



(B8) 



RA 



Here we have introduced the source angle 4> defined by C = sin(0/2). The matrices V{x) and V~^{x') operate in the 
retarded-advanced space only and hence commute with any disorder operators placed between the Green functions. 
As a result, factors V and drop from expressions for any closed diagrams. The matrices A in the above equation 
allow us to decompose the Green function into a direct product of the two operators acting in the RA space and in 
the sublattice space. The matrix A commutes only with the potential disorder and must be retained as a part of the 
Green function in the general case. 

The regularized Green function arising in the calculation of diagrams with point-like impurities takes especially 
simple form when V{x) and l/^^(a;) are singled out, 



Greg (a;) = lim 

a;'— fx 



U[x,x,y)+^^ (a;-x')2 + 2;2 



/ 1 



V 



\ 

'<Jx— \ 



sin^y 



V-\x). 



(B9) 



It is worth noting that the above result does not depend on the order of taking the limits. 



Appendix C: Resonances in p-wave scattering 

In order to include higher scattering resonances in our consideration, we take into account small deviations of the 
regularized Green function endpoints from the center of impurity. Expansion up to the first order in these deviations 
allows for the resonances in the p scattering channel. Let us define the following matrix containing the regularized 
Green function and its various derivatives: 



5(l'o) 



lim 

r->-ro 
r'->-ro 



1 dx' 
dx dxdx' 

dy dydx' 



dxdy> 
dydy 



Go(r,r')-g(r-r') 



(CI) 



Here and below we use the matrix notation with brackets rather than parentheses to distinguish the newly introduced 
matrix from the internal structure of the Green function in the RA and sublattice spaces. The [1, 1] component of 
Greg coincides with the regularized Green function Greg from Eq. (|B9I) . 
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The impurity is described by its T matrix which is an operator with the fohowing kernel: 
r(r,r') = C/(r),5(2)(r - r') + U{v)g{v - v')U[v') 

oo „ 

+ Y, rf'ri . . . d?Vp U{v)g{v - ri)U{v,)g{r, - r^) . . . U(vp)g[vp - v')U{v'). (C2) 
p=i"' 

In the main text we have focused on the s-wave resonances and characterized the impurity by its integrated T matrix. 
In the simplest case of circular potential impurity the integrated T matrix is nothing but the scattering length i in 
the s channel. Here we retain more information about the shape of the impurity and introduce the following integral 
of the T matrix: 



T 



1 x' -Xq y' - yo 

x-xo [x - xo){x' - xo) {x - xo){y' ~ Vo) 

y-yo {y - yQ){x' - xq) {y - yo){y' - yo) 



T(r,r'). 



(C3) 



Here and below we assume r ~ {x, y} and r' = {x' , y'}. The [1, 1] component of T is the integrated T matrix used in 
the main text. 

With the above-defined matrices 
the form 



Jrcg and T, we can express the single-impurity correction to the free energy in 
ST = In dot [1 - TGrcg(ro)] . (C4) 



This expression includes the product of [1, 1] components of Gicg and T, reproducing the result obtained in the main 
text, along with the next terms of the Taylor expansion of G{r, r') in the vicinity of Tq. 

Let us now calculate the matrix Greg- Taking derivatives of G(r,r') we reveal the following internal structure of 



Greg (a:) 



1 

ay 

(7rr 



Greg (a:) Grog (2;) 

^ice{x) Greg (a;) 





'l 





0" 







ay 





(C5) 



Here the element Greg is given by Eq. (jB9[) and the other elements are 

a2 ^2 



Gl,ix)^—Vix) 



Glo,ix)^^V{x 



/ 30" - it' 



cos 
''^ sin^ ^ 



-la,, 



\ 



-la 
^2 



^TTsin^ 



TT sin ^ 
3</.2 - cos ^ 

Cz -r—^ h i(T. 



7r0 COS 



67r2 



a„:„3 ■KX_ I 

T ' 



sm 



TTc/) COS ^ 



\ sin ^ 



7r2 + 02 



Y-\x\ 



V-\x). 



(C6) 



(C7) 



sin^f; 2sin^ / 



Taking the advantage of the representation (|C5P we can reduce the size of the matrix in Eq. (jC4[) by introducing 



1 

ay ax 



T 



1 
ay 
ax 



1 

ay{x - xo) ^ax{y- yo) 



T(r,r') 



ay{x' -xo) + ax{y' - yo) 



t 



Now the correction to the free energy takes the simplified form 

,5J" = lndctil-T 



Grog(a;o) GJ'(,g(xo) 

Greg(a;o) G,*eg(a;o), 



(C8) 



(C9) 



In order to obtain the final result from Eq. (|C9p we have to establish a relation between the matrix T and the 
solution of the scattering problem with a given impurity profile. Let us place the impurity at the origin Tq = and 
consider the scattering state with the following asymptotic behavior far from the impurity: 



ikr 



kr > 1. 



(CIO) 
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This expression contains an incoming plane wave with the momentum k and an outgoing spherical wave with the 
direction-dependent scattering amplitude fk- We use the notation tp^ i^r) for the angle between the vector k (r) and 
the direction of x axis. 

The scattering amplitude is directly related to the matrix clement of the T matrix between the incoming and 
outgoing waves 



fk{f,(p') = ^ j cPvcPv'il e r(r,r') ^^,y^ exp [ifc(a;'cos(/7' + 2/'sin(p' - xcos.^ - ysinv?)]. 



(Cll) 



The integrals entering the matrix T can be expressed through the fc — > limit of the scattering amplitude and its 
derivatives in the following way: 



lim 

fc-i-O 



2-rr 



dipd(p' ( 1 e""^' 



27r2 



,itp i{tf>-tf>') 



1 



-4i su\(2(p') 



Ok 



4isin(2(p)^ 8 sin(2(p) sin(2(p'^ ^' 



(C12) 



Once the scattering amplitude for a given impurity is known, one can use Eqs. (jC12[) and (IC9[) to find the correction 
due to such an impurity to the full counting statistics, including s- and p-wave scattering channels. Below we apply 
this general result to the circular potential impurity. 

In the case of axially symmetric impurity, the scattering amplitude depends on the difference f — ^' and Eq. (|C12p 
simplifies to 



t: VO e"^ 



1 

4cos(2^)^ 



fk{^p)- 



(C13) 



We can get rid of the ip integral and express the result in terms of the scattering phases 6m where m is a half-integer 
number — z component of the total momentum (angular -|- spin) of the Dirac electron. The scattering amplitude is 
the following Fourier series with the coefficients determined by the scattering phases: 



2i5„ 



gi(m-l/2)^^ 



(C14) 



In the case of potential impurity preserving the time-reversal invariance, the scattering phases obey the symmetry 
relation 5m = S-m- At low energies they decay according to Sm ~ 0(fc2|™l). Using these properties, we substitute 
Eq. (|C14P into Eq. (|C13|) and obtain 



'£ 
h 



^ = -41im 

fc-s-o k ' 



-16 lim 



33/2 



(C15) 



Note that, contrary to the parameter ix has a dimension of the third power of length. 

If the potential impurity has the rectangular profile with the height mq and radius a, the Dirac equation is easily 
solved yielding the following scattering parameters: 



27ra 



Ji(uoa) 



27ra 



3 J2(uoa) 



Jo(uoa) Ji(uoa) 
Now we evaluate the determinant in Eq. (|C9|) and obtain the result 

2 



(C16) 



1 



16L2 7r2 



Jl TTXq 



256i6 



128L4 



- 1 



97r2 



1 



_,2 -Kxa 



4cos2 

„■ 6 ££0 " 



sm 



+ cot 



2 ^ 

L 



^o{ee^\. (C17) 



Here we have neglected the term ^ since it is always much smaller than the term ^ U\. In Fig. [5] we compare 
the conductance calculated from Eq. (jC17l) as a function of u^a with the results of numerical simulations [25j. 
Exactly at the p resonance, when — >■ cxd, correction to the generating function becomes 



(5 J" 2 In 



62 2 



97r2 



4 cos2 ^ 



(CIS) 
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FIG. 5: Correction to the conductance from a circular impurity placed at Xo = L/2 as a function of uoa. Dotted line: numerical 
simulations 25] with the parameters a/L — 0.2, W/L — 6; solid line: analytic result including s and p resonances. 

From this expression we find the correction to the conductance and average it with respect to the impurity position 



8e2 f dx d^ST 



/^sin^x(lB + siny e^^ 



(C19) 



This value is bigger than the conductance correction at the s-wave resonance, Eq. ([T}. 



Appendix D: Two-impurity correction 



The second-order correction to the free energy involves the diagrams with two impurities shown in Fig. [S] We sum 
up these diagrams assuming the two impurities are located at the points {x,yQ} and {x',yQ + y} and obtain 



llndetjl- [l-TG,eg(x)] ^TG{x,x' ;y)[l - TG^.^ix')] Vg(x', x; -j/)}. 



In the resonant case T — > oo, this correction to the free energy takes the form 

1 



■ hi dot 



1 - G;,lix)Gix,x';y)G~lix')Gix',x;-y) 



With the Green functions calculated above, we evaluate the determinant and obtain the result 
(2) , (j}"^ - (j}'^ {a\ + + b\ + b'^_ + cl + d) + (a+a_ + b+b- + c+c_) 

(02_4)(02_„2_) 



a± 



sin7r(u ± v) 



b+ = 



TT 



X + X 

2L 



sin7r(f ± w) 
X — x' 

V 



sin7r(?x; ± u) ' 



2L 



w 



2L' 



(Dl) 
(D2) 

(D3) 
(D4) 

(D5) 



This expression is to be averaged with respect to the impurity positions. To calculate corrections to the conductance 
and noise, we expand the free energy in small 4> up to the forth order and integrate each term over x, x' , and y 
numerically. Together with the zeroth- and the first-order terms in riimp, this yields 



7" = 



4^ 



+ 4:nin,pWL In (vr + x/tt^ - cj)^) + nf^^^W / dx dx' / dy S^^'>T 

Jo J -oo 



= const — - 



W 



47r ^ V 



4 

87r4 



- nLpi^ ( 0.1290 2 



0.08823 



24 



(D6) 
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FIG. 6: Feynman diagrams for one- and two-impurity corrections to the free energy. 



The conductance and the Fano factor up to the second-order corrections are 



h 86^ 



=0 

1 2 d*T/d^%=o 1 12 r^^n^ao-, 2 ,4 



4 16 



1.032 nf^pL^ 



(D7) 
(D8) 



